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^ ' Abstract 

^ , Using the analysis of ChPT to two loops, we perform an error analysis of the 

Oh' threshold and low energy parameters, based on the uncertainties for the one loop 



low energy parameters and the resonance saturation mechanism. Different sets of 
t"^ . one loop low energy constants have been considered. Thus, the predictive power 

of the effective field theory, is quantified on the basis of the present experimental 
uncertainties. 
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Qv ! Chiral perturbation theory (ChPT) is the effective pion field theory where the expan- 

sion parameter is m? / {\f^ fY , with m the pion mass and / the weak pion decay constant. 
P^l At some finite order an increasing number of undetermined parameters is generated which 

p ■ can, so far, only be fixed by experimental data. Besides the theoretical uncertainty intro- 

^ • duced by the finite number of terms considered in the expansion, there is an additional 

experimentally induced uncertainty. The low energy parameters inherit such an uncer- 
tainty, thus limiting in practice the predictive power of the effective field theory; if the 
r> I errors at a given order, say (9(m^"/(v^47r/)^"') are larger than the contributions of the 

c3 ■ next order 0(m^"'"'"^/(\/47r/)^""'"^), the calculation of the latter becomes useless, unless 

more accurate experiments are performed. 

To make our points quantitative, we consider tttt scattering as the prototype reaction to 
study the role played by chiral symmetry breaking and the validity of the chiral expansion 
within QCD at low energies. The scattering lengths and effective ranges for the isospin 
I = 0,1,2 have been computed in SU{2) ChPT at tree level |I|, one loop 0| and two 
loops 0, 3, in a power series expansion 

au = 47 + a}^°°P + a?i°°P + ■ ■ ■ (1) 

hjj = b]Y + bYr^ + h^r" + ■■■ . (2) 

At tree level, the number of parameters involved are two: the pion weak decay constant 
/tt = 93.2 MeV and the pion mass m = 139.6 MeV. For the purposes of the present 
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discussion the experimental error bars in these parameters can effectively be taken to be 
zero. At the one loop level one has four, in principle undetermined, parameters Ji, I2, I3, h- 
At the two loop level the whole SU{2) vrvr amplitude information can be gathered into six 
coefficients, which in Ref. 0, || have been called 61, . . . , feg- Ideally, these parameters could 
be extracted from a direct low energy analysis of vrvr scattering experiments. The data 
are, however, too poor in the low energy region and several methods have been devised. 
Motivated by the success of the resonance saturation hypothesis at scales /i ~ 0.5 — 1 
GeV, at the one loop level P, the authors of Ref. ||^ suggest to estimate the bulk of 
the two loop corrections through resonance saturation. This approach requires the values 
for the one loop /'s deduced from some clean source. In Ref. 0] two different sets of I 
parameters have been considered: i) by taking the /i and I2 parameters obtained from K14 
form factors using a dispersive one loop calculation for three flavors (Set I in ||^), and 
ii) fixing /i and I2 to the values deduced from D— wave vrvr scattering lengths, at order 
0{p^) (Set II in [^). Unfortunately, no error analysis has been undertaken in Ref. [^. 
On the other hand, the estimates based on Kn analysis, can now be improved thanks to a 
recent study of this process at two loop accuracy ^. Other possibility is to extract the one 
loop low energy constants, together with their errors, from Roy sum rules, saturated by 
the high energy behavior of vrvr scattering [§]. In conjunction with resonance saturation, 
this procedure also yields an error analysis of threshold parameters . 

An error analysis from the point of view of the predictive power of ChPT is missing, 
and it is the main subject of the present work. We assume that primary quantities, 
coming from a x^~fit or direct experimental measurements, are Gauss distributed (with 
or without correlations), as we learn from elementary statistics. We propagate errors 
by means of Monte Carlo simulations, keeping always statistical correlations between 
all parameters entering in a given derived quantity. We use 10^ samples, and to quote 
errors we always use a 68% confidence level around the central value. Since the out- 
coming threshold parameters distributions are not gaussians, we take in this way into 
account possible skewness in the distributions. By using the Monte Carlo method we 
avoid summing errors in quadratures, which would be incorrect for statistically correlated 
quantities, and we do not have to use any complicated covariance formula. 

Let a*^""-* and Aa*^"-' be the n-loop central value and error of an observable. Thus, to be 
predictive and convergent at the ra— loop level one ought to have the relation, 

Aa(") « a("+^) « a("\ (3) 

In order to establish these necessary requirements from a quantitative point of view we 
will adopt the currently accepted scheme of assuming that resonance saturation yields the 
bulk of the two loop contribution at scales 0.5 — 1.0 GeV. This induces a scale ambiguity 
in the two-loop part of the 6— parameters, which we implement by providing the scale 
with an error, 

/i = 750 ± 250 MeV . (4) 

Moreover, since one does not expect resonance saturation to be exact, we provide these 
parameters with 100% uncertainty and take, in the notation of Ref. Q, the values 

n = -0.6 X 10"^ X (1±1) ra = 1.3 X 10"^ X (1 ± 1) rg = -1.7 x 10"^ x (1 ± 1) 
r4 = -1.0 X 10"^ X (1 ± 1) rs = 1.1 x 10"^ x (1 ± 1) rg = 0.3 x 10~^ x (1 ± 1) (5) 



at the scale /z given in Eq. (^). In addition to these numbers, one needs estimates of the 
one loop low energy parameters ^1,2,3,4- Through all this work we fix 

[3 = 2.9 ± 2.4 , U = 4:A ± 0.3 , (6) 

as determined from the study of SU{3) breaking effects and the scalar form factor in 
Refs. and respectively. There is less consensus regarding the values of Ji and I2, 
therefore we consider three possibilities. We give below the central values, errors and 
linear correlation coefficients for /i and I2 in each case, in order to characterize the full 
statistical distributions, deduced from different primary distributions, which will be used 
later to calculate threshold and low energy parameters. 

• Set I Ki4 form factors: Ref. allows to determine h and I2 from a two loop analysis 
of Ki4 decays. The analysis is done in terms of the scale dependent SU{3) low energy 
constants L[ at the p mass scale, rrip = 770 MeV, and central values for /i and I2 are 
given. The relation between the SU{2) and SU{3) low energy constants is known 
at one loop [jlO[ and thus the determination of the h and I2 parameters in suffers 
from two loop SU{3) uncertainties. From the remarks of that work a strong anti- 
correlation between li and I2 can be inferred, although the corresponding correlation 
matrix is not given. We will take here two extreme view points, total de-correlation 
and total anti-correlation. As we will see, the totally anti-correlated case seems 
to induce smaller errors in the scattering lengths as compared to the totally de- 
correlated case. We will call these two choices Set la and Set lb respectively. We 
will also examine here a partial anti-correlation scenario (Set Ic) which models more 
accurately the limited statistical information provided in Ref. 0. More details are 
given in the Appendix. 

For the three sets, we take L[ = 0.52 ±0.23, L^ = 0.72 ±0.24 and L^ = -2.69 ±0.99 
(main fit in Table 1 of Ref. [0]). For Set la, we assume that the three statistical 
distributions are de-correlated, whereas for Set lb [Set Ic] assume these are corre- 
lated as specified in Eq. (|^) [ Eq. (|13D with ri3 = r23 = — ri2 = —0.85]. By means 
of a Monte Carlo simulation we generate the statistical distributions for the /i and 
I2 parameters, and thus we get 



h = 0.3 ±2.1, 


fa = 4.77 ±0.45, 


r{hj2) = 


0, 


Set la 


(7) 


h = 0.3 ±1.0, 


fa = 4.77 ±0.45, 


r{liM) = 


-1, 


Set lb 


(8) 


h = 0.3 ±1.2, 


[2 = 4.77 ±0.45, 


r{kM) = 


-0.69, 


Set Ic 


(9) 



being r the linear correlation coefficient. The errors quoted above are clearly low 
bounds, because they do not account for any systematic effects, in particular those 
induced by the (9 (p®)— corrections to the relations between the two- and three- 
flavor low energy constants^. Estimates for the systematic errors are not given in 
Ref. either. Notice that correlations among L\2z *io '^o^ affect to the error in I2 
(see footnote 4). 



'^Note that in the second entry of Ref. 0| there is an inconsistency. In Eq. (6.24) an error of ±1.0 for T2 
is quoted. That error is not compatible with having I2 — 1927r^L2 ~ (1 + lii("^/i'//^^) + ^ ln(m^//i^))/8, as 
deduced from Eq. (6.23) of that reference, and an error for Lj of ±0.24 x 10^'^. The previous formula gives 



Set II JJ— waves: The method of Ref. allows to define another parameter set. 
If the D— wave vrvr scattering lengths for isospin 7 = and / = 2, ao2 and 022 
respectively, are fixed at the two loop level one gets parameter Set II, 

^1 = -0.8 ±4.8, /2 = 4.5 ±1.1, r(Ji,l2) = -0.75 , Set II (10) 

To obtain /i and I2 and their errors, we have propagated the errors in /i, 002, 022, 
I3, U and the resonance parameters of Eq. (^) in the formula for the D— wave 
scattering length given in Ref. [§]. This procedure generates a correlation between 
/i and I2 which has to be taken into account when calculating errors in quantities 
depending on the previous parameters. Note that our central values for Zi and I2 
numbers are not exactly the ones quoted in Ref. |Q since they take /i = 1 GeV. 
Furthermore, in Ref. @] no error estimates are quoted for /i and /2- If we had taken 
H = 1.00 ± 0.25 GeV, we would have obtained h = —1.5 ± 5.8 and I2 = 4.5 ± 1.1, in 
agreement with their quoted central value. 

Set III Roy sum rules: Finally, using the method of Roy sum rules, Ref. [^ , another 
parameter set has been obtained. We call it Set III, 

[1 = -0.9 ±1.2, [2 = 4.34 ±0.25, r(ri, [2) = -0.22 , Set III (11) 

Both central values and errors quoted above do not agree with those given in Eq. (21) 
of Ref. P]. The difference on the central values is due to the fact that in this 
latter reference all resonance parameters are set to zero at 1 GeV (rj(l GeV) = 
± 2 X 10"'^). To obtain the errors we have propagated the errors in the resonance 
saturation scale /i, Ai, A2 (parameters given in Eq. (2) of Ref. [§), Is, [4 and the 
resonance parameters of Eq. (|^) when Eqs.(6-7) of Ref. ^ are inverted. This, 
again generates a correlation between li and I2 as in the previous casesQ. The Ai, A2 
parameters suffer from sizeable systematic uncertainties. Those have been estimated 
in the second entry of Ref. and turned out to be of comparable magnitude, 
when not bigger, than the statistical fluctuations, and as a consequence statistical 
correlations get washed out. This justifles the use of de-correlated distributions for 
these two parameters both in Ref. M and in the present work. 



an error for I2 of ±1927r^ x 0.24 x lO^'^ — ±0.45 and it is obviously independent of the possible correlations 
between the SU{3) L'[,L2 and L^^parameters, because it only involves L^. This is in agreement with 
our results. The problem in Ref. m is that the error definition for the 5't/(2)— parameters, "projections 
on the relevant variable of the 68% confidence level domain" , is not consistent with that adopted for the 
^[/(S)— parameters. For these latter ones, the standard and traditionally accepted x^— errors are given. 
The method of the projections would lead to significantly different errors for L^ 2 3, as can be appreciated 
from the figures shown in Ref. m. In anycase, this projection method is not standard. The standard and 
complete procedure consists of giving the full correlation matrix deduced from the x^— fit. We ignore to 
what extent the errors on derived quantities given in Ref. ||7| are affected by this inconsistency. 

^Errors differ from those quoted in Ref. Q, because of the different treatment considered here (error 
on the resonance saturation scale, different choice, both for central values and errors, for the resonance 
parameters but also because there is a numerical mistake in that work, which we correct here. The error 
analysis of Eqs. (16) and (17) of that work yields h = -0.37 ± 0.95 ± 1.27 and h = 4.17 ± 0.19 ± 0.33, 
which errors do not agree with those quoted in Eq. (20) of that reference. 
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Figure 1: Two dimensional distributions of the low energy parameters li,h generated via Monte 
Carlo simulation with N — 10"* samples. The contours and the line made of squares for the Set lb 
represent the 68% confidence limit, i.e., inside them the 68% of the samples are enclosed . Upper left 
panel, parameter Sets la and lb corresponding to a two loop Ki4^ form factor analysis 0, with totally 
de-correlated and totally anti-correlated ^i and I2 parameters respectively. Upper right panel, parameter 
Set Ic corresponding to a two loop Ki4^ form factor analysis 0, with linear correlation r = —0.69. Middle 
left panel, parameter Set II corresponding to D— waves in tttt scattering (r = —0.75). Middle right panel, 
parameter Set III corresponding to Roy sum rules in tttt scattering (r = —0.22). Finally, in the bottom 
panel the 68% confidence limit contours are superposed. 



The determination of the one loop parameters, their corresponding uncertainties as well 
as their mutual correlation can be best seen by plotting the resulting two-dimensional 
Monte Carlo generated distributions for all considered sets. This is done in Fig. |l]. For all 
sets we consider a sufficiently large sample, A^ = 10^, and have checked numerical stability 
by doubling the sample size. We also plot the 68% confidence limit contours, i.e., inside 
these boundaries 68% of distribution points are enclosed. For a better comparison we have 
also superposed (bottom panel) the 68% confidence limit contours to check for statistical 
compatibility. As could already be anticipated, the main difference between parameter 
Sets la, Ic and lb (corresponding to K14 decays) is that the errors decrease as the anti- 
correlation increases. The parameter Set II provides the largest uncertainties, and large 
overlap with the remaining sets is encountered. At the same time, it is clear that there 
is a trend to discrepancy between parameter Sets Ic and III, since very little overlap is 
observed between both distributions. We should also mention that the level of uncertainty 
regarding both the K14 and the Roy sum rules determinations is comparable, although 
the latter method seems to be slightly more accurate. 

In Ref. 0] the calculation of the two loop contributions to the threshold parameters 
was undertaken and explicit expressions for scattering lengths and ranges were written, 
in terms of the low energy constants, &i,2, 3,4,5,6- In this note we complete their numerical 
calculation by providing their numbers with the inherited error-bars and also exploring 
different parameter sets. This is shown in Table |1] (for the fejparameters) and in Table ^ 
(for the threshold parameters, aij,hjj). 

For all sets, we have generated statistical distributions both in the one loop and the 
two loop corrections to the scattering lengths and effective ranges by propagating the 
statistical distributions in the low energy parameters /i^2,3,4, the resonance parameters 
^1,2,3,4,5,6 and the scale /x disccused above. The uncertainty in the scale and in the resonance 
contributions, affect to the error of the two loop corrections only. Errors in the two loop 
contributions are correlated with those in the one loop contributions. Therefore, the error 
of the sum cannot be obtained by simply adding the errors in quadrature of one and two 
loop contributions. In any case, for the two loop calculation to be numerically meaningful 
these one loop uncertainties have to be significantly smaller than the corrections due to 
the two loop calculation. At the same time the two loop correction has to be significantly 
smaller than the one loop correction, for a convergent expansion. 

As we see in Table |1], all predictions for the hi parameters are compatible within 
errors. Moreover, as one might have inferred from the values and distributions of the 
parameters h and T2 the increasing accuracy of the determinations correspond to Sets II, 
la, Ic, lb and III, in that order^. Besides, we would like to point out that for the Sets 
Ia,Ib and Ic the statistical fiuctuations in the one loop parameters induce errors on the 
hi parameters as important, when no more as in the de-correlated case of Set la, as the 
scale uncertainty quoted in Ref. 0]. For the parameter Set II the central values of the h 
parameters are slightly different than those in Ref. [Q because the central li and I2 values 
are also different. The errors for the parameter Set II are much larger than those found 



^ These two latter parameter sets have similar accuracy. 
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-0.27 + 0.65 


-0.27 i °:1I 


10^ ■ 64 
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Table 1: Low energy scale independent two loop parameters and their uncertainties due to the error bars 
in the I's one loop parameters and the uncertainties in both the scale and the resonant part of the two 
loop contribution. Set la corresponds to K14 two loop calculation: li — 0.3 ± 2.1 , I2 = 4.77 ± 0.45 , I3 — 
2.9 ± 2.4 , I4 = 4.4 ± 0.3, with li and I2 uncorrelated. Set lb corresponds to K14 two loop calculation: 
h = 0.3 ± 1.0 , h = 4.77 ± 0.46 , [3 = 2.9 ± 2._4 , k = 4.4 ± 0^3, with h and r2_totally anti-correlated. Set 
Ic corresponds to K^ two loop calculation: h = 0.3 ± 1.2 J2 = 4.77 ± 0.44 , I3 = 2.9 ± 2.4 , U = 4.4 ± 0.3, 
with a correlation coefhcient r{li,l2) — —0.69. Set II corresponds to a two loop 13— wave tttt scattering 
lengths calculation (see main text): li = -0.8 ± 4.8 , h = 4.5 ± 1.1 , I3 = 2.9 ± 2.4 , I4 = 4.4 ± 0.3, with a 
correlation coefficient r(li,l2) — —0.75. Set III corresponds to a Roy equation sum rule analysis of tttt 
scattering at two loops: h = -0.9 ± 1.2 , ^2 = 4.34 ± 0.25 , [3 = 2.9 ± 2.4 , U = 4.4 ± 0.3, with a correlation 
coefficient r(li, I2) — —0.22. The scale is /i = 750 ± 250 MeV under the resonance saturation hypothesis, 
for which a 100% error is assumed. 



for parameter Sets I and III. The size of the errors for Set II is comparable with the ones 
for 63, 64, 65, be found in Ref. IQjQ 

As we see in Table 0, results for 020, 620 and 022 with parameter Set la are such that 
the one loop errors are larger than the central values of the two loop contribution. The 
inclusion of statistical correlations, Sets lb and Ic, increases the accuracy of the predic- 
tions, though it still happens that errors on the one loop contribution make irrelevant 
the two loop one. The situation worsens dramatically, for parameter Set II where we see 
that in most considered cases predictive power is lost beyond one loop, with the exception 
of the S"— wave scattering lengths aoo and 092. Set III, obtained from the Roy-sum rule 
analysis of Ref, [Q, turns out to be as predictive as Sets lb and Ic, obtained from the 
two-loop improved Ki, analysis of Ref. 0. In any case is also true that for 020, 620 and 
022 within Set III the one loop uncertainties are larger than the two loop contribution. 
The results of the table are compatible, for all parameter sets with the experimental anal- 
ysis of TTTT scattering data []rT| but produce in general much better errors despite of the 

^ Regarding the parameter Set III, one should say that if one would calculate the 6's from Eq. (48) 
of Ref. H], with the numerical input of Eq.(2) of Ref. M, one would have obtained 

10 5i = -0.68+0.09 , 10 62 = +0.64 ±0.09 , 10^63 = -0.35 ± 0.24 
10^64 = +0.47 ±0.03 , 10^65 = 0.13 ±0.06 , lOHg = 0.10 ± 0.01 

&3.4,5.6 are in agreement with Eq.(49) of Ref. j^. We sec that this is not the same as evaluating the &'s 
in the spirit of resonance saturation (Eqs. D.1,D.2 and D.3 of second entry in Ref. [Q) from previously 
computed values of Zi and I2, through their relation to Ai and A2 parameters, as suggested in Ref |8| and 
adopted in this work. In any case, the values given above and those corresponding to Set III in Table |l| 
are compatible within statistical fluctuations. 
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Figure 2: Left panel: contours corresponding to a 68% confidence limits for the two dimensional 
distributions of the S'-wavc scattering lengths aoo and 020- Parameter sets are defined in the main text. 
The linear correlation coefficients, r(aoo,a2o), for Sets la, lb, Ic,II and III are 0.85,0.38,0.58,0.86 and 
0.56 respectively. Right panel: The same as before with the additional inclusion of the universal bands 
and the 70% confidence limit contour (crosses) deduced from a numerical analysis of the Roy equations 
in Ref. fO. 



problems discussed above on the relevance of the two loop contributions. The early work 
on Ki4 form factors using a dispersive one loop calculation for three flavors , predicted 
a isoscalar D— wave three standard deviations above the value extracted from the experi- 
ment [||, the upgrade of this calculation to two loop accuracy in Ref. have contributed 
to considerably improve such a discrepancy. 

In Fig. 1^ we exhibit the Monte Carlo propagated two dimensional distributions of 
the S'-wavc scattering lengths, for isospin channels / = and 1 = 2, for all parameter 
sets and in terms of the 68% confidence limit contours. It is interesting to compare this 
distributions with the preliminary results of Ref. [|r^, based on the detailed numerical 
analysis of Roy equations in vrvr scattering. Such an analysis implies the existence of a 
universal band, outside which no solution to Roy equations exists, as a mathematical 
consequence of analyticity, crossing, unitarity and isospin conservation. In addition to 
this, the physical requirement of describing the available experimental data, provides a 
70% confidence limit contour [0]. As can be seen in the figure, with all parameter sets 
considered in this work, the two-loop ChPT analysis are not only compatible with the 70% 
confidence limit contour of Ref. [Q but have significant smaller statistical fluctuations. 



In summary, effective field theories like Chiral Perturbation Theory have predictive 
power, but it is not unlimited because of three reasons: 1) truncation of the expansion, 2) 
proliferation of undetermined constants at any order and 3) experimental uncertainties. 
Thus, experimental data prove crucial to determine the, increasing with the order, un- 
known constants and their errors, which propagate in a correlated way to higher orders in 
the expansion possibly undermining the "convergence" of the expansion. This situation 
also appears in fundamental theories like QED or QCD, but it is in fact worse in ChPT 



because the number of unknown parameters in these theories does not increase with the 
order of the expansion. This why error analysis is so important. We have exemphfied 
our points in the calculation of the threshold parameters for vrvr scattering up to two 
loops. The general picture provided by ChPT is rather satisfactory, in the sense that 
the accuracy of the predictions is much bigger than the available data Q. Nevertheless 
we find cases where, within the present experimental accuracy, the errors in the one loop 
contribution are larger than the central values of the two loop contribution. Conclusions 
regarding the loss of predictive power, can only be reinforced if systematic uncertainties 
induced by the the 0(p^)— corrections to the relations between the two- and three-flavor 
low energy constants are included in the results of Ref . [0 . The effect of these systematic 
corrections will be twofold: a general increase of the error fluctuations and a decrease of 
the correlations. Both effects will contribute to increase the errors on all derived threshold 
and low energy parameters presented in this work and deduced from i^i4— decays. 

ChPT defines a whole family of effective theories, but obviously the most interesting 
choice exactly corresponds to low energy QCD. In QCD with two flavors and neglecting 
isospin breaking there appear only two independent parameters, Aqcd and the quark mass 
iriQ. The quantities /^r, rn, the Vs, the 6's and higher order parameters must be functions 
of them. This dependence introduces correlations among all low energy parameters which, 
if taken into account, would influence the present error analysis, and presumably might 
yield to more moderate errors in terms of Aqcd and rriQ and their errors. Assuming more 
statistically independent parameters than QCD suggests is, so far, another manifestation 
of the inability to undertake a quantitative and microscopic derivation of ChPT as an 
effective low energy theory of QCD on the one hand, but on the other hand reassures ChPT 
as a convenient tool to deal with non-perturbative phenomena in strong interactions. 



Appendix: Modeling statistical correlations in K14 decays. 



In Ref. is stated that L3 is strongly anti-correlated with LI and L2, though central 
values and errors for these parameters are given (L^ = 0.52 ± 0.23, L2 = 0.72 ± 0.24 and 
L3 = —2.69 ± 0.99, main fit in Table 1 of Ref. [0), however the correlation matrix is not 
provided. Besides, new parameters Xi = Lg — 2L^ ^L^, X2 = L^ and X3 = (Lg — 2L^)/L3 
are introduced, it is quoted the value X3 = 0.12l^^ and it is also specified that X3 is little 
correlated with Xi and X2. We will show below how statistical modeling and the fragmen- 
tary information given in Ref. [0 can be used to reconstruct some relevant information 
on the correlations between li and /2- To be most objective, we also explore two other 
extreme cases which correspond to total de-correlation and total anti-correlation. This 
said, we explore three different scenarios in this work: 

• A total de-correlated picture {r{L'[,Ll) = r{L2,L^) = r{L[,L2) = 0, being r the 
linear correlation coefficient). It leads to X3 = O.I2I20) ^"(-^15-^3) = —0.14 and 
r(X2,X3) = -0.09. 



^Right plot of Fig.|| constitutes a clear example. 



A total anti-correlation scenario {r{Ui, L^) = r{L2,Ll) = —1) which can be imple- 
mented by simply assuming a linear relation 

m - {L\)) /a, = - {L; - {LD) /a, = {U, - {U,)) /a,, (12) 

That is the necessary and sufficient condition for total anti-correlation. Besides, one 
also has r^L^^U^) = +1- This simple model leads to X3 = 0.12^6 and r(Xi,X3) = 
r(X2, X3) = 0.06. Thus, the total anti-correlation scenario, though simple, provides 
an acceptable description, and in any case much more precise than when correlations 
between L^, Lg and L3 are neglected, of the findings of Ref. |0]. 

A partial anti-correlation scenario which can be implemented if one assumes that 
L\ , L2 ^iid L3 are gaussian distributed^ according to 

P(L;.L^.« = (detPI)^ .eri^c-.) (13, 

(27r) 2 0-10-20-3 

where P is the join density probability distribution of the three random variables, 
cTj is the Ul error, L is a column matrix with elements (L[ — fii)/<7i, being /ij the 
central value of the variable L[, and finally the inverse of the symmetric matrix C 
is given by 

/ 1 ri2 ri3 \ 
C-' = ri2 1 r23 (14) 

V ^^13 r23 1 / 

and therefore it is determined by the linear correlation coefficients r^j = r{L^,Uj). 
Our aim is to improve the total anti-correlated scenario presented in the item above 
but keeping the model still as simple as possible. Thus we have explored situ- 
ations where ri3 = r23 = — T12 = r, for which one can analytically diagonalize 
the C— matrix in Eq. ([T3|). We find for r = —0.85, still a strong anti-correlation, 
X3 = 0.12t?o, r(Xi,X3) = -0.02 and r(X2,X3) = -0.01. Thus, this model for the 
correlations reproduces pretty much the results given in the K14 analysis of Ref. [^ . 
This can also be appreciated in Fig. |^ where we show the resulting two dimensional 
distributions of the SU{3)—\ow energy parameters L[,Ll and L2,Ll. These dis- 
tributions compare reasonably well with those given in Fig.2 of the ffist entry of 
Ref. [^ . This makes us more confident on the validity of the simple statistical model 
used here. 



The correlations between L\, Lg and L3 are transported through the equations relating 
SU{2) and SU{3) low energy constants. For the total anti-correlated scenario leads to a 
total anti-correlation for li and I2, it is to say (Ji — (li))/ai = —{h — {l2))/cr2 and therefore 
^^(^1,^2) = —1- For the total de-correlated case the correlation coefficient r(li,l2) is zero, 
whereas for the partial correlation scenario we get r{li,l2) = —0.69. 

^This is totally justified, because in Ref. m these variables have been determined from a x^— fit. 



10 




CJ 




3 "-3 

Figure 3: Two dimensional distributions of the SU{3)— low energy parameters L\,L'^ and L^^L^- 
generated via Monte Carlo simulation with N — l{f samples. The contours represent the 68% confidence 
limit, i.e., inside them the 68% of the samples are enclosed. 
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o/j; bij 


Set 


(tree) 


+ (lloop) 


±(21oop) 


total 


experiment 


aoom 


la 


0.156 


0.043 ±0.005 


0.015 ±0.004 


0.214 ±0.008 


0.26 ±0.05 




lb 


U 


" ± 0.002 


" ±0.003 


" ± 0.004 


u 




Ic 


u 


" ± 0.003 


" ±0.003 


" ± 0.005 


u 




II 


u 


0.040 ±0.007 


0.013 ±0.004 


0.209 ±0.010 


u 




III 


u 


0.039 ±0.003 


0.012 ±0.002 


0.208 ±0.005 


u 


boomP 


la 


0.179 


0.070 ±0.015 


0.025 ±0.011 


0.273 ±0.024 


0.25 ±0.03 




lb 


u 


" ± 0.003 


u + 0.007 
- 0.008 


u + 0.008 
- 0.009 


u 




Ic 


u 


" ± 0.006 


±0.008 


u + 0.011 
- 0.012 


u 




II 


u 


0.059 ±0.024 


0.019 ±0.013 


r. 9C-7 + 0.032 
V.ZOI _ 0.036 


u 




III 


u 


0.058 ±0.008 


0.018 ±0.006 


0.254 ±0.010 


u 


10 ■ aiim^ 


la 


0.297 


0.055 ±0.012 


0.021 t °:°°1 


n Q7Q + 0.011 
U.-J'J - 0.012 


0.38 ±0.02 




lb 


u 


" ± 0.009 


u + 0.002 
- 0.003 


" ±0.011 


u 




Ic 


u 


" ± 0.009 


u + 0.002 
- 0.003 


" ±0.011 


u 




II 


u 


0.059 ±0.033 


0.018 1 °:°°I 


n Q7P; + 0.035 

yj.oio _ 0.031 


u 




III 


u 


0.059 ±0.007 


0.018 ±0.002 


0.374 ±0.007 


u 


10 ■ bum^ 


la 





0.030 ±0.012 


n nop: + 0.006 
U.UZO _ g 008 


0.055 t n?? 


— 




lb 
Ic 




u + 0.008 
- 0.009 

" ± 0.009 


u + 0.002 

- 0.004 
u + 0.003 

- 0.005 


u + 0.009 

- 0.011 
u + 0.009 

- 0.011 


: 




II 


u 


0.034 ±0.033 


0.020 t S 


0.054 ±0.029 


— 




III 


u 


0.034 ±0.007 


0.019 1 ',Zl 


n ncrq + 0.005 
U.UOO „ 0.006 


— 


10 ■ 020"^ 


la 


-0.446 


0.021 ±0.021 


0.002 ±0.004 


-0.423 ± 0.024 


-0.28 ±0.12 




lb 


ii 


" ± 0.008 


" ± 0.004 


" ±0.010 


a. 




Ic 


u 


" ±0.011 


" ± 0.004 


" ±0.013 


u 




II 


u 


0.007 ±0.028 


0.000 ±0.004 


-0.439 t ',Z 


u 




III 


cc 


0.004 ±0.013 


-0.001 t S 


-0.442 ±0.015 


C4 


10 ■ 620"^^ 


la 


-0.892 


0.129 ±0.043 


0.003 1 HI? 


-0.760 ± 0.043 


-0.82 ±0.08 




lb 
Ic 




" ±0.010 

u + 0.018 
- 0.019 


u + 0.012 
- 0.010 
u + 0.013 
- 0.011 


u + 0.016 
- 0.014 
u + 0.023 
- 0.021 


a 




II 


u 


0.095 ± 0.045 


0.003 + °:°li 


-0.794 ±0.039 


u 




III 


u 


0.088 ±0.022 


0.003 ±0.011 


-0.801 ±0.021 


u 


10^ ■ 002"^^ 


la 





0.143 ±0.031 


0.058 t °:°1^ 


0.202 ±0.042 


0.17 ±0.03 




lb 
Ic 




" ± 0.009 
" ±0.015 


u + 0.016 
- 0.019 

u + 0.016 
- 0.019 


u + 0.021 

- 0.025 
u + 0.024 

- 0.028 






II 


u 


0.117 ±0.028 


0.053 ±0.022 


0.170 ±0.030 


u 




III 


u 


0.111 ±0.016 


0.051 1 'ofil 


0.163 ±0.020 


u 


10^ • a22m^ 


la 





0.278 ±0.241 


-0.035 1 H? 


0.24 ±0.17 


0.13 ±0.30 




lb 
Ic 




" ± 0.062 
" ±0.108 


u + 0.057 
- 0.039 

u + 0.064 
- 0.044 


u + 0.07 

- 0.05 
cc + 0.10 

- 0.08 






II 


u 


0.119 ±0.460 


0.011 t till 


0.13 ±0.30 


u 




III 


cc 


0.102 ±0.130 


0.015 t Ifi 


0.12 ±0.11 


cc 



Table 2: Threshold tttt scattering parameters and their uncertainties in units of m due to the error bars 
in the I's one loop parameters and the uncertainties in both the scale and the resonant part of the two 
loop contribution. Sets la, lb, Ic, II and III are defined in table 1 and main text. Errors are not added 
in quadrature due to statistical correlations. Experimental values are from Ref. pj . 
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